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Abstract.  The  response  of  a  general  circulation  model  to  a  change  in  its  treatment  of 
cloud  solar  forcing  is  investigated.  Radiation  field  data  from  the  forecast  model  of  the 
Navy  Operational  Global  Atmospheric  Prediction  System  for  five  Julys  (1979-1983)  are 
presented  in  an  investigation  of  the  effect  of  a  change  from  grid  cell  averaged  clouds  to 
maximally  overlapping  clouds  'n  the  model's  solar  radiation  scheme.  The  model  results 
are  compared  with  Nimbus  7  Earth  Radiation  Budget  top  of  the  atmosphere  (TOA) 
solar  and  longwave  irradiances  and  with  derived  surface  solar  irradiance  data. 

Although  the  maximal  overlap  scheme  performs  considerably  better  than  the  grid  cell 
averaging  scheme  (reducing  maximum  deficiencies  in  TOA  and  surface  solar  irradiance 
by  over  100  W  m“'),  significant  errors  remain.  The  simulated  correlation  between  TOA 
net  solar  and  longwave  irradiance  improves  at  low  latitudes  in  the  northern 
hemisphere,  with  little  change  at  higher  latitudes.  This  improved  correlation  is 
consistent  with  the  greater  consistency  between  the  treatments  of  solar  and  longwave 
cloud  radiative  forcing  brought  to  the  model  by  the  new  solar  radiation  scheme.  The 
change  in  the  radiation  treatment  is  shown  to  have  the  greatest  direct  effect  on  solar 
radiation  over  convective  regions,  a  consequence  of  the  scarcity  of  optically  thick 
clouds  produced  by  the  model’s  cloud  parameterization  in  other  regions.  The  model 
responds  with  an  increase  in  convective  activity  over  land  and  an  increase  in  the  flux 
of  moisture  from  sea  to  land  Planetary  cooling  over  the  oceans  increases  because  of  a 
decrease  in  cloud  cover.  From  mid  to  high  latitudes  in  the  northern  hemisphere,  there 
are  scattered  regions  of  increased  cloud  water  content  associated  with  increased 
tropospheric  temperatures.  Over  land  the  model  response  in  terms  of  TOA 
downwelling  solar  irradiance  tends  to  counter  the  increase  in  solar  irradiance  caused  by 
the  model  change  in  all  latitudinal  zones  in  the  northern  hemisphere.  This  response  is 
caused  primarily  by  changes  in  the  cloud  fields,  which  thus  act  as  a  negative  feedback 
following  the  change  in  cloud  solar  forcing.  The  significance  of  this  response  is 
examined  with  respect  to  the  perturbation  in  solar  irradiance  represented  by  the  model 
change.  An  estimate  of  this  perturbation  is  obtained  by  taking  the  difference  in  solar 
irradiance  diagnosed  by  the  two  cloud  solar  forcing  treatments  for  simulations 
employing  the  grid  cell  averaging  scheme.  The  response  is  significantly  greater  in 
magnitude  in  the  tropics  than  at  midlatitudes,  both  in  an  absolute  sense  and  as  a 
percentage  of  this  perturbation.  Because  TOA  longwave  irradiance  exhibits  a  positive 
response  in  the  tropics,  and  a  negative  response  at  midlatitudes,  however,  the 
percentage  response  in  net  TOA  downwelling  irradiance  is  actually  greater  in 
magnitude  at  midlatitudes.  In  a  number  of  regions  the  cloud  feedback  is  very  large, 
showing  the  importance  for  cloud  field  prediction  of  improvements  in  the  treatment  of 
cloud  solar  forcing.  Such  cloud  feedback  also  explains  the  small  improvement  seen 
here  in  the  prediction  of  TOA  solar  irradiance  in  certain  regions.  Increases  in  surface 
sensible  heating  and  longwave  cooling  are  generally  considerably  less  than  increases  in 
surface  latent  heating,  though  a  notable  exception  occurs  in  arid  central  Asia.  A  large 
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ground  temperature  increase  in  that  region  is  strongly  correlated  at  low  levels  with  the 
atmospheric  temperature  increase  observed  at  midlatitudes  in  the  northern  hemisphere. 


1 .  Introduction 

The  importance  of  clouds  in  modulating  the  radiative 
heating  of  the  Earth-atmosphere  system  has  stimulated  a 
substantial  amount  of  research  into  cloud  radiative  effects  in 
recent  years.  Research  on  global-scale  cloud  nidiative  ef¬ 
fects  has  included  both  observational  studies  (e.g.,  London, 
19.57;  Budyko.  1969;  Cess.  1976;  Oliriiif;  and  Clapp.  1980; 
Ramanathan  el  al..  1989;  Ardaniiy  el  a!..  1991;  Smith  and 
Yonder  Haar.  1991;  Hartmann  et  al..  1992;  Laszlo  and 
Pinker,  I993|  and  modeling  studies  [e.g..  Cess  and  Potter, 
1987;  Slinyo  and  Slinpo.  1988;  Harshvardhan  et  al..  1989; 
Cess  et  al..  1990;  Morcrette,  1990;  Slin.i’o  and  Slingo,  1991; 
Potter  et  al..  1992).  To  gain  insight  into  how  model  accuracy 
can  be  improved,  a  growing  number  of  studies  have  com¬ 
pared  model  radiation  fields  with  satellite  observations  or 
with  satellite-derived  d:oa  (e.g  .  G"lpyr.  e*  u/.,  19^7;  Slingo 
et  al.,  1987;  Buriez  et  al.,  1988;  Kiehl  and  Ramanathan. 
1990;  Bony  et  al..  1992;  Chertoek  and  Sad,  1^^' 

Because  of  the  complex  nature  of  clouds  number  of 
uncertainties  are  generally  associated  with  modeling  studies 
of  cloud  radiative  properties.  Variations  in  the  microphysical 
structure  of  clouds  can  affect  cloud  radiative  properties 
tremendously  [e.g..  Wiseomhe  et  al..  1984;  Ackerman  and 
Stephens.  1987;  Lion  et  al..  1991 1.  as  can  variations  in  cloud 
cover,  cloud  macrophysical  shape,  and  the  spatial  distribu¬ 
tion  of  water  in  clouds  [e.g..  McKee  and  Cox.  1974;  Aida. 
1976;  Davies.  1978;  Weinman  and  Harshvardhan,  1982; 
Harshvardhan.  1982;  Welch  and  Wielicki.  1984;  Kite.  19871. 
Fhe  complex  nature  of  clouds  is  particularly  problematic  for 
attempts  to  simulate  their  effect  on  the  large-scale  circula¬ 
tion.  General  circulation  models  (GCMs)  are  constrained  at 
present  to  using  only  very  crude  treatments  of  clouds  and 
their  radiative  interactions.  In  such  models,  cloud  properties 
arc  often  parameterized  in  terms  of  large-scale  variables, 
though  explicit  schemes  that  carry  a  prognostic  variable  for 
cloud  water  have  also  been  developed  [e.g..  Smith.  1990).  In 
the  diagnostic  schemes,  cloud  amount  is  assumed  to  depend 
on  such  quantities  as  large-scale  relative  humidity,  convec¬ 
tive  precipitation,  vertical  potential  temperature  gradient, 
and  vertical  velocity  [e.g..  Slingo.  1987).  Radiative  transfer 
in  cloudy  grid  cells  is  generally  treated  using  calculations 
that  are  based,  at  least  to  some  degree,  on  the  (plane- 
parallel)  assumption  that  clouds,  and  hence  radiation  fields, 
do  not  vary  at  all  in  the  horizontal. 

Our  work  examines  the  effects  of  a  change  from  the  use  of 
grid  cell  averaged  clouds  to  maximally  overlapping  clouds  in 
the  solar  radiation  scheme  of  the  spectral  foreca.st  model  of 
the  Navy  Operational  Global  Atmospheric  Prediction  Sys¬ 
tem  (NOG APS).  The  assumption  of  maximal  cloud  overlap 
is  employed  in  the  model's  longwave  radiation  scheme 
\Harshvardhan  et  al..  1987]  but  has  not  previously  been  used 
in  the  NOCiAPS  solar  radiation  scheme,  where  it  is  poten¬ 
tially  considerably  more  computer  intensive  than  the  grid 
cell  averaging  treatment  currently  used.  Further  discussion 
of  the  NOCiAPS  GCM  is  given  in  section  2  along  ^ 
(llscussion  uf  iiic  simuiaiiv^;-,.,  tliat  were  used  for  the  present 
work.  In  the  first  part  of  our  study  the  model  results  are 
compared  with  Nimbus  7  Earth  Radiation  Budget  (ERB) 


observations  in  order  to  assess  the  net  effect  of  the  change  in 
cloud  solar  forcing  with  respect  to  observed  model  errors. 
This  comparison  is  presented  in  section  4.  following  a 
discussion  in  section  3  of  the  radiation  data  that  are  used. 
Both  the  direct  effect  of  the  change  in  the  treatment  of  cloud 
solar  forcing  and  the  model  response  to  that  change  in  terms 
of  radiation  fields  are  examined  in  section  5.  In  section  6  the 
response  of  model  cloud  fields  is  examined.  Changes  in 
precipitation  and  ground  and  atmospheric  temperatures  are 
discussed  in  section  7. 

2.  Model  and  Simulation  Details 

Simulation  data  from  version  3.3  of  the  NOGAPS  forecast 
model  (NOGAPS  3.3),  taken  from  a  10-year  simulation 
beginning  January  1979.  were  used  for  this  work.  This 
simulation  was  carried  out  at  the  Naval  Research  Labora¬ 
tory  (NRL)  as  part  of  the  Program  for  Climate  Model 
Diagnosis  and  Intercomparison  (PCMDI).  Data  from  simu¬ 
lations  for  the  five  consecutive  Julys  beginning  in  1979, 
carried  out  with  a  maximal  overlap  version  of  the  NOGAPS 
model,  described  below,  were  also  used.  Fields  from  the 
NOGAPS  3.3  PCMDI  simulations  were  used  to  initialize  the 
model  for  these  simulations,  and  a  1 -month  spin-up  period 
was  used  for  each  simulation. 

Version  3.3  of  the  NOGAPS  forecast  model  is  a  spectral 
model  with  18  vertical  levels.  This  model  was  implemented 
for  the  PCMDI  simulations  with  a  triangular  truncation  at 
wavenumber  47.  which  corresponds  roughly  to  a  2.5°  trans¬ 
form  grid.  The  radiation  scheme  is  that  of  Davies  [1982]  and 
Harshvardhan  et  al.  [1987],  the  solar  portion  of  which  is  an 
extension  of  the  scheme  developed  by  Lads  and  Hansen 
[1974].  The  only  gases  assumed  to  De  absorbers  of  solar 
radiation  are  water  vapor  and  ozone,  and  the  effects  of 
aerosols  are  ignored.  The  cloud  scheme  is  similar  to  the 
diagnostic  parameterization  of  Slingo  [1987].  The  prediction 
of  convective  clouds  differs  from  that  of  Slingo,  however,  in 
that  convective  precipitation  at  the  surface  is  used  to  diag¬ 
nose  convective  cloud  cover,  rather  than  the  precipitation  at 
cloud  base  level.  More  detailed  descriptions  of  the  model  are 
given  by  Hogan  and  Rosmond  [1991],  who  present  a  discus¬ 
sion  of  NOGAPS  3.2,  and  by  Hogan  and  Brody  [1993],  who 
list  changes  incorporated  into  version  3.3  of  NOGAPS. 
These  changes  include  the  adoption  of  the  assumption  of 
maximally  overlapping  clouds  in  the  model's  treatment  of 
longwave  radiation  and  the  change  to  a  more  realistic  value 
of  single  scatter  albedo  for  radiative  transfer  of  solar  radia¬ 
tion  in  clouds.  The  parameterization  of  cloud  optical  thick¬ 
ness  that  was  used  for  the  PCMDI  simulations  is  also  a  new 
feature  of  the  model.  The  present  parameterization  was 
developed  by  Tim  Hogan  at  NRL  and  is  described  in  the 
appendix.  Although  this  parameterization  appears  to  pro¬ 
duce  optical  thicknesses  that  are  too  low  in  the  case  of 
convective  clouds,  and  will  tend  to  reduce  to  some  extent 
differences  b'-' v  een  the  two  radiation  schemes  examined 
here,  this  fact  should  not  significantly  affect  the  results  of  our 
present  study. 

Where  partial  cloudiness  occurs,  the  treatment  of  long¬ 
wave  radiation  in  NOGAPS  3.3  follows  the  maximal  overlap 
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treatment  described  by  Harshvardhan  et  al.  [1987],  There  is 
no  provision  in  the  Harshvardhan  et  al.  radiation  scheme, 
however,  for  a  similar  treatment  of  solar  radiation.  In  the 
solar  radiation  scheme,  partially  cloudy  grid  cells  are  treated 
using  a  grid  cell  averaging  treatment  in  which  the  optical 
thickness  of  clouds  is  obtained  by  averaging  cloud  optical 
thickness  over  the  entire  grid  cell.  This  was  not  much  of  a 
problem  in  the  University  of  California,  Los  Angeles/Goddard 
Laboratory  for  Atmospheres  (UCLAi/GLA  GCM  for  which 
the  Harshvardhan  et  al.  scheme  was  designed,  because 
partial  cloudiness  was  allowed  only  for  optically  thin  bound¬ 
ary  layer  clouds  [Harshvardhan  et  at.,  1989].  In  NOG.V  'S 
3.3,  however,  the  inconsistency  between  the  treatment  »f 
partial  cloudiness  used  for  solar  radiation  and  that  used  for 
longwave  radiation  is  potentially  of  considerable  conse¬ 
quence.  We  found  that  this  inconsistency  can  be  avoided  in 
an  economical  manner  by  repeated  application  of  the  Harsh¬ 
vardhan  et  al.  scheme.  In  this  treatment,  solar  irradiances 
are  summed  for  each  different  vertical  profile  of  cloudiness 
in  a  given  grid  column  consistent  with  the  assumption  of 
maximal  cloud  overlap  (see  Figure  1).  The  revised  model  in 
which  this  maximal  overlap  treatment  is  used  for  solar 
radiation,  denoted  here  as  NOGAPS  3.3MO,  requires  only 
about  109?-  more  CPU  time  on  a  CRAY  Y-MP  than  does 
NOGAPS  3.3. 

3.  Radiation  Data 

The  model  simulation  results  were  compared  with  July 
mean  top  of  the  atmosphere  (TOA)  net  downwelling  solar 
and  outgoing  longwave  irradiances  for  the  years  1979-1983 
derived  from  wide-angle  field-of-view  (WFOV)  Nimbus  7 
ERB  Experiment  observations  [Jacohowitz  et  at.,  1984;  Kyle 
et  at..  1985]  and  with  corresponding  July  mean  surface  solar 
irradiances  over  the  oceans  derived  from  Nimbus  7  ERB 
planetary  albedo  data  by  Chertock  et  al.  [1992]  using  the 
model  of  Chertock  [1989]  and  Frouin  and  Chertock  [1992]. 
These  data  are  based  on  observations  from  a  single  Sun- 
synchronous  satellite,  which  means  that  modeled  diurnal 
cycles  of  albedo  and  outgoing  longwave  radiation  had  to  be 
used  to  produce  the  monthly  mean  data.  It  has  been  shown, 
however,  that  the  ERB  TOA  data  compare  favorably  with 
Earth  Radiation  Budget  Experiment  (ERBE)  data  derived 
from  observations  from  a  trio  of  sun-synchronous  satellites: 
ERBS.  NOAA  9.  and  NOAA  10  [Kyle  et  at..  1990].  The 
RMS  differences  between  WFOV  ERB  and  ERBE  data 
interpolated  to  the  higher-resolution  2.5°  grid  of  the  ERBE 
scanner  data  for  April  and  July  1985  are  approximately  4  W 
m  ■  for  longwave  irradiances  and  8-9  W  m  ’  for  solar 
irradiances  at  the  top  of  the  atmosphere  [Kyle  et  at.,  1990]. 
The  RMS  differences  between  the  WFOV  ERB  data  and  the 
ERBE  scanner  data  for  the  same  months  and  gnu  aic 
approximately  10  W  m  ’  for  longwave  irradiances  and  14-15 
W  m  '  for  solar  irradiances.  Although  the  level  taken  to  be 
the  TOA  in  processing  the  ERB  data  was  somewhat  arbi¬ 
trarily  taken  to  be  15  km  above  sea  level,  it  has  been  shown 
that  raising  this  level  to  30  km  above  se.a  level  p:  or'y 

a  fraction  of  a  percent  difference  in  calculated  WFOV 
longwave  irradiances  [Kyle  et  at..  1990].  The  surface  solar 
irradiance  data  used  here  have  been  estimated  to  be  accurate 
to  within  10-20  W  m  '  for  monthly  averages  at  the  9°  grid 
resolution  of  the  present  study  [Chertock  et  at..  1992]. 

The  availability  of  data  near  the  poles  varied.  TOA  net 


Figure  1.  Vertical  cross  section  through  a  model  grid  col¬ 
umn.  The  maximal  overlap  assumption  is  implemented  for 
the  hypothetical  cloudiness  profile  shown  here  by  applying 
the  model’s  radiation  scheme  separately  to  each  of  the  four 
indicated  vertical  columns. 


solar  irradiances  (obtained  from  albedo  and  insolation  data) 
were  available  between  72°S  and  72°N  for  the  five  Julys.  The 
same  was  generally  true  of  outgoing  longwave  irradiances, 
except  for  July  1980.  for  which  data  were  not  avail.able 
poleward  of  54°N.  The  surface  solar  irradiance  data  were 
available  between  63°S  and  63°N. 

4.  Error  in  Model  Radiation  Fields 

4.1.  Solar  Radiation  Comparisons 

The  mean  TOA  net  downwelling  solar  irradiance  observed 
from  ERB  data  for  the  five  consecutive  Julys  is  shown  in 
Figure  2a.  The  difference  between  the  corresponding  values 
obtained  with  the  models  and  the  observed  values  is  shown 
in  Figure  2b  for  NOGAPS  3.3  and  in  Figure  2c  for  NOGAPS 
3.3MO.  Such  differences  will  be  referred  to  hereinafter  as 
errors,  though  they  are  not  necessarily  true  errors  because  of 
inaccuracies  in  the  ERB  and  derived  surface  solar  irradiance 
data.  Shading  is  used  here,  and  in  all  our  difference  maps,  to 
indicate  where  the  Student’s  t  test  assigns  a  significance  to 
the  plotted  mean  differences,  scaled  with  the  pooled  monthly 
variability  of  these  means,  at  the  95%  level.  The  results  in 
Figure  2  show  a  striking  improvement  with  the  maximal 
overlap  scheme,  though  substantial  errors  still  exist.  Over 
the  Pacific  Ocean,  for  example,  the  maximum  deficiency  of 
TOA  net  solar  irradiance  obtained  with  the  NOGAPS  model 
is  reduced  from  172  W  m  ■  to  49  W  m~"  with  the  new 
scheme.  Both  model  versions  produce  large  errors  over  the 
Indian  Ocean,  though  NOGAPS  3.3MO  does  somewhat 
better.  The  peak  deficiency  in  that  region  is  reduced  from 
152  W  m  '  to  91  W  m  ““  with  the  maximal  overlap  scheme. 

The  results  in  Figure  2  show  a  large  amount  of  regional 
variability  in  the  differences  betwc'-n  the  grid  cell  averaging 
and  the  maximal  overlap  treatments.  In  order  to  better 
understand  the  observed  differences,  we  show  in  Figure  3a 
the  mean  surface  convective  rainfall  rate  from  the  NOGAPS 
3.7  siiuulaiioiis  foi  iiiv  fi.c  Julys.  A  comparison  ot  the 
convective  rainfall  shown  in  Figure  3a  with  the  TOA  solar 
irradiance  error  shown  in  Figure  2b  indicates  a  strong 
correlation  between  simulated  convective  rainfall  and  defi¬ 
ciencies  in  TOA  net  solar  irradiance  predicted  by  the  cloud 
grid  cell  averaging  scheme  of  NOGAPS  3.3.  In  Figure  3b. 
which  shows  the  differences  between  the  convective  rainfall 
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a)  ERB 


b)  NOGAPS  3.3  -  ERB 


Figure  2.  (a)  July  mean  TOA  net  solar  irradiance  (in  watts  per  square  meter)  from  ERB  measurements, 

(bl  difference  between  results  from  NOG.APS  3..^  and  the  data  shown  in  Figure  2a,  and  (c)  difference 
between  results  from  NOGAPS  3.3MO  and  the  data  shown  in  Figure  2a.  Shaded  regions  are  those  for 
which  the  plotted  differences  are  statistically  significant  at  the  95'7r  level  according  to  the  results  of  a 
Student's  t  test.  The  contoui  interval  is  20  W  m"’.  and  the  zero  contour  has  been  om'tted.  Negative 
contours  are  indicated  by  dashed  lines. 


rate'^  predicted  by  the  two  model  versions,  we  see  that 
althougti  there  are  some  significant  fractional  increases  in 
convective  rainfall  over  land,  particularly  over  Africa  and 
North  and  South  America,  differences  ir  convective  activity 
do  not  explain  the  large  differences  between  Figures  2a  and 
2b.  The  increased  convection  obtained  "''h  NOC.'.PS 
3.3,N10  u,cr  land  should,  in  tael,  tend  to  increase  the 
deficiencies  in  TOA  solar  irradiance.  but  this  is  not  ob¬ 
served.  The  correlation  between  the  error  in  Figure  2b  and 
the  convective  rainfall  simulated  by  NOGAPS  3.3  appears  to 
stem,  rather,  from  the  optically  thick  clouds  parameterized 
in  regions  of  deep  convection  by  the  cloud  parameterization 
scheme  used  in  the  NOG.APS  model.  For  water  clouds. 


optical  thickness  is  assumed  to  be  proportional  to  the 
saturation  specific  humidity,  and  it  is  the  relatively  low  base 
levels  of  convective  clouds  that  are  most  responsible  for 
their  significant  optical  thickness  in  the  NOGAPS  simula¬ 
tions.  The  poor  performance  of  the  grid  cell  averaging 
-chCiTiC  foi  k-onvective  djuJ:.  most  likely  stcm,>  froui  .he 
well-known  fact  that  for  optically  thick  plane-parallel  clouds, 
albedos  asymptotically  approach  a  limiting  value  (e.g.. 
Chandrasekhar.  I960,  pp.  89-91:  Arking  and  Childs.  1985: 
Stephens  and  Tsay.  1990].  When  one  averages  cloud  optical 
thickness  over  model  grid  cells  for  optically  thick  clouds,  the 
effect  of  the  decrease  in  the  albedo  of  the  cloud-covered 
region  can  be  smuh  in  comparison  with  the  net  increase  in 
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Figure  3.  (a)  July  mean  surface  convective  precipitation  rate  (in  millimeters  per  day)  predicted  by 
NOGAPS  3.3  and  (b)  the  difference  between  the  corresponding  results  obtained  with  NOGAPS  3.3MO  and 
the  results  obtained  with  NOGAPS  3.3.  The  contour  interval  is  2  mm  d  “ ' .  Shading  is  defined  as  in  Figure 
') 


the  grid-scale  albedo  that  results  from  having  a  greater  area 
covered  by  cloud.  An  illustration  of  this  fact  is  shown  in 
Figure  4,  where  the  albedos  obtained  using  the  two  cloud 
solar  forcing  schemes  considered  here  for  a  cylindrical  cloud 


0  2  .4  8  .8  1 

Cloud  Fraction 

Figure  4.  Albedo  plotted  against  cloud  fraction  for  a  cylin¬ 
drical  cloud  of  an  optical  thickness  of  328  that  extends  from 
1.1  km  above  the  surface  to  a  height  of  6.5  km. 


of  optical  thickness  equal  to  328  that  extends  from  1.1  km 
above  the  surface  to  a  height  of  6.5  km  are  plotted  as  a 
function  of  cloud  fraction.  The  results  shown  for  the  grid  cell 
averaging  scheme  are  the  same  as  would  be  obtained  for  a 
cloud  fraction  of  unity  and  a  cloud  optical  thickness  that  is 
varied  linearly  from  zero  to  328.  The  grid  cell  averaging 
scheme  thus  produces  an  asymptotic  dependence  of  albedo 
on  cloud  fraction  that  seems  considerably  less  reasonable 
than  the  linear  relationship  obtained  with  the  maximal  over¬ 
lap  scheme. 

Although  the  maximal  overlap  scheme  in  NOGAPS 
3.3MO  yields  significantly  improved  results  over  regions  of 
deep  convection,  significant  errors  remain,  both  in  convec¬ 
tive  regions  and  elsewhere.  We  note  here  the  model's 
tendency  to  predict  excessively  large  values  of  TOA  net 
solar  irradiance  over  regions  where  extensive  marine  strati¬ 
form  cloud  coverage  should  occur.  Where  low-level  stratus 
or  stratocumulus  clouds  generally  predominate  in  nature,  a 
substantial  portion  of  the  model  error  apparently  results 
from  the  fact  that  the  model  significantly  underestimates  the 
occurrence  of  such  clouds  \Hogan  and  Brody.  1993].  For 
July  1985,  for  example,  NOGAPS  predicts  a  cloud  fraction 
of  0.1  (defined  here  to  be  the  greatest  cloud  fraction  at  any 
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a)  ERB 


b)  NOGAPS  3.3  -  ERB 


c)  NOGAPS  3.3M0  -  ERB 


Figure  5.  (a)  July  mean  surface  solar  irradiance  (in  watts  per  square  meter)  derived  from  ERB 

measurements,  (b)  difterence  between  the  results  from  NOGAPS  3.3  and  the  data  show  n  in  Figure  5a,  and 
(c)  difference  between  the  results  from  NOGAPS  3.3MO  and  the  data  shown  in  Figure  5a.  Shading  and 
contours  are  defined  as  in  Figure  2. 


model  grid  level)  for  a  5“  latitude  (I2.5'’S  to  I7.5°S)  by  7.5° 
longitude  (80°W  to  87.5°W)  region  off  the  west  coast  of  South 
.America  known  to  be  dominated  by  stratus  clouds.  For  the 
same  region  and  month,  Ockcrt-BfU  and  Hartmann  [1992] 
recently  reported  a  low-level  (below  680  mbar)  cloud  frac¬ 
tion  of  0.6.  Similar  model  deficiencies  have  been  reported  in 
some  related  model  validation  studies  [Bitriez  el  a!..  1988; 
Vesperini  et  al..  1991 ;  Soden,  I992|.  The  failure  of  NOG.APS 
in  this  regard  stems  from  its  requirement  for  an  unrealisti¬ 
cally  strong  low-level  inversion  for  stratus  clouds  to  occur 
(T.  Rosmond.  personal  communication.  1993).  a  problem 
which  has  since  been  corrected. 

The  5-year  July  mean  surface  solar  irradiance  data  are 
presented  in  Figure  5a.  The  plotted  values  appear  to  be 
rather  highly  correlated  with  the  values  of  TO.A  net  down- 
welling  solar  irradiance  plotted  in  Figure  2.  as  similarly- 
noted  by  Laszlo  and  Pinker  [19931.  The  differences  between 
the  model  predictions  and  the  ERB-derived  data  are  shown 


in  Figures  5b  and  5c.  Over  regions  of  considerable  convec¬ 
tive  activity  the  maximal  overlap  scheme  yields  a  substantial 
improvement  in  predicted  values  of  surface  solar  irradiance 
similar  to  that  shown  in  Figure  2  for  TOA  solar  irradiance. 
Peak  deficiencies  of  194  W  m  '  and  96  W  m  "  are  obtained 
for  NOGAPS  3.3  and  NOGAPS  3.3MO.  respectively;  the 
former  value  occurred  in  the  central  Pacific  Ocean,  and  the 
latter  value  occurred  in  the  Indian  Ocean.  The  deficiency  of 
194  W  m  ■  in  the  central  Pacific  obtained  with  NOGAPS  3.3 
is  reduced  to  53  W  m  "in  the  NOGAPS  3,3MO  simulations. 

Comparisons  of  the  model-generated  zonal  means  of  TOA 
net  solar  irradiance  with  ERB  data  are  shown  in  Figure  6. 
Over  many  latitudinal  zones,  regional  errors  in  NOG.APS  3.3 
tend  to  compensate  to  produce  reasonably  accurate  zonal 
means.  .Such  is  not  the  case  over  the  northern  hemisphere 
tropics,  however,  where  large  deficiencies  of  net  solar 
irradiance  exist.  The  revised  model  does  not  produce  the 
sizable  errors  in  tropical  zonal  means  that  characterize 
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Figure  6.  Zonal  July  mean  TOA  net  solar  irradiance.  Re¬ 
sults  from  NOGAPS  3.3  (plus  signs),  and  NOGAPS  3.3MO 
(solid  circles)  are  shown  with  ERB  data  (open  circles). 


NOGAPS  3.3.  but  it  docs  exhibit  a  general  tendency  to 
overpredict  the  amount  of  net  downwelling  solar  radiation, 
particularly  for  midlatitudinal  and  high-latitudinal  zones  of 
the  northern  hemisphere.  Some  of  liiis  error  results  from  the 
failure  of  the  model  to  predict  enough  low'-level  marine 
stratus  clouds,  as  indicated  :n  the  discussion  above.  The  fact 
that  NOGAPS  3.3MO  tenos  to  err  somewhat  more  on  the 
high  side  than  does  NOGAPS  3.3  may  well  result  from  an 
inadequacy  of  the  maximal  overlap  assumption  for  certain 
regions,  though  this  conclusion  cannot  be  positively  deduced 
from  our  results. 

4.2.  Longwave  Radiation  Comparisons 

A  plot  of  mean  TOA  outgoing  longwave  irradiance  ob¬ 
tained  from  the  ERB  Experiment  for  the  five  Julys  is  shown 
in  Figure  7,  together  with  plots  of  the  difference  between  the 
model-derived  values  and  the  observations.  The  results 
indicate  considerable  errors  in  the  longwave  irradiance  field 


a)  ERB 
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c)  NOGAPS  3.3M0  -  ERB 


Figure  7.  (a)  July  mean  TOA  outgoing  longwave  irradiance  (in  watts  per  square  meter)  from  ERB 

measurements,  (b)  difference  between  the  results  from  NOGAPS  3.3  and  the  data  shown  in  Figure  7a.  and 
(c)  difference  between  the  results  from  NOGAPS  3.3MO  and  the  data  shown  in  Figure  7a.  Shading  and 
contours  are  defined  as  in  Figure  2. 
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Figure  8.  July  mean  TOA  net  soliir  irradiance  tor  9"  x  9" 
regions  extending  from  0’  to  9'N  plotted  against  July  mean 
TOA  outgoing  longwave  irradiance:  (a)  results  from 
NOC’iAPS  J.3  (plus  signs)  and  from  F.RB  measurements 
(circles)  and  (b)  results  from  NOC/AP.S  3.3MO  (plus  signs) 
and  from  ERB  measurements  (circles). 


in  several  regions  for  both  model  versions.  The  models  err 
on  the  low  side  over  the  Indian  Ocean,  for  example,  by  as 
much  as  72  W  m  *  for  the  NOGAPS  .T.^  simulations  and  74 
W  m  ’  for  the  NOCiAPS  .^._3MO  simulations.  For  several 
regions,  including  the  Indian  Ocean  and  the  central  Pacific 
Ocean,  a  comparison  with  Figure  .3  shows  that  such  low 
values  of  predicted  outgoing  longwave  irradiance  coincide 
with  large  amounts  of  simulated  convective  rainfall.  That 
this  correspondence  i .  not  universal  is  seen  in  the  results 
obtained  for  .Australia.  Indonesia,  and  the  tropical  eastern 
Pacific  Ocean.  Off  the  west  coasts  of  North  America  and  the 
northern  portion  of  South  .America,  the  model  predicts  too 
much  outgoing  longwave  radiation.  T  his  observation  is  con- 
si' tent  with  the  general  failure  of  the  model  to  predict 
realistic  amounts  of  low-level  marine  stratiform  clouds. 

4.3.  Relationship  Between  Solar  and  Longwave 
Radiation  Fields 

To  accurately  simulate  the  net  radiative  forcing  by  clouds, 
it  is  reasonable  to  assume  that  a  model  must  provide 
accurate  forecasts  id' the  correlation  between  longwave  and 
solar  radiation  fields.  In  order  to  assess  the  effects  of  the 
model  change  studied  here  in  this  respect,  scatterplots  of 
rOA  net  solar  irradiance  versus  TO.A  outgoing  longwave 


irradiance  for  NOG.APS  .3. .3  and  NOG.APS  .3..3.V10  are  pre¬ 
sented  in  Figure  8  along  with  the  corresponding  F.RB  data. 
The  data  show  n  are  for  9’  x  9^^  regions  that  extend  from  the 
equator  to  9'N.  We  see  a  considerable  improvement  in  going 
from  the  clo’id  grid  cell  averaging  scheme  of  NOG.APS  3.3  to 
the  maximal  overlap  scheme  of  NOG.APS  .3.3MO,  The 
apparent  limitation  of  the  cloud  parameterization  is  mani¬ 
fested  in  both  plots,  however,  in  some  very  low  values  of 
outgoing  longwave  irradiance,  the  occurrence  of  which  has 
already  been  noted  in  Figure  7.  The  results  shown  here 
support  the  conclusion  that  the  substantial  improvement  in 
the  solar  irradiance  field  over  convective  regions  obtained 
with  the  maximal  overlap  scheme  (as  seen  in  Figures  2  and  .‘3) 
does  not  result  primarily  from  a  compensation  by  that 
scheme  for  errors  in  the  parameterized  cloud  fields.  Corre¬ 
sponding  plots  for  midlatitudes  to  high  latitudes  show  con¬ 
siderably  less  difference  between  the  two  model  versions,  a 
reflection  of  the  fact  that  the  models  differ  the  most  in  the 
regions  of  greatest  convective  activity. 

5.  Response  of  Model  Radiation  Fields 

5.1.  Direct  Effect  on  the  Absorption  of  Solar  Radiation 

The  comparisons  presented  in  section  4  show  that  the 
change  in  cloud  solar  forcing  examined  here  has  a  consider¬ 
able  effect  on  NOGAPS  TOA  and  surface  radiation  fields. 
The  observed  changes  may  be  attributable  not  only  to  the 
direct  effect  of  the  change  in  the  treatment  of  cloud  solar 
forcing  but  also  to  the  response  of  the  model  cloud  fields 
(and  to  a  lesser  extent  the  moisture  field)  to  this  change.  In 
order  to  better  understand  the  significance  of  the  observed 
model  response,  we  first  address  the  question  of  the  extent 
of  the  direct  effect  represented  by  the  model  change.  .An 
estimate  :if  the  magnitude  of  this  effect  was  obtained  by 
repeating  the  NOGAPS  .3. .3  simulations  for  the  five  Julys 
using  the  same  cloud  grid  cell  averaging  scheme  that  was 
used  for  the  original  simulations  and  outputting  diagnostic 
solar  radiation  fields  calculated  using  the  maximal  overlap 
scheme. 

The  model  response  in  the  present  simulations  derives 
from  the  effect  of  the  model  change  on  the  amount  of  solar 
radiation  absorbed  by  the  atmosphere  and.  for  land  regions, 
by  the  surface.  (Because  sea  surface  temperatures  are  set  to 
observed  values,  the  large  change  in  surface  solar  irradiance 
over  the  oceans  seen  in  the  difference  plots  in  Figure  has 
no  effect  on  the  simulations.)  The  relative  contribution  of 
each  of  these  effects  is  illustrated  in  Figure  9.  w  hich  show  s 
the  difference  w  ith  respect  to  solar  radiation  absorbed  by  the 
surface  (Figure  9a)  and  by  the  atmosphere  (Figure  9b) 
between  values  obtained  from  solar  irradiance  fields  diag¬ 
nosed  from  the  NOG.APS  .3.3  July  simulations  using  the 
maximal  overlap  scheme  and  fields  based  on  the  cloud  grid 
cell  averaging  scheme  for  the  same  simulations.  The  amount 
of  solar  radiation  found  to  he  absorbed  by  the  surface  using 
the  maximal  overlap  scheme  is  generally  considerably 
greater  than  that  obtained  using  the  grid  cell  averaging 
scheme.  This  is  tnie  p;irticularly  over  tropical  convective 
regions,  as  is  similarly  observed  in  FTgure  .x.  The  amount  of 
solar  radiation  found  to  be  absorbed  by  the  atmosphere  using 
the  two  treatments  qf  cloud  solar  forcing  differs  by  rather 
modest  amounts,  with  the  grid  cell  averaging  scheme  yield¬ 
ing  more  absorption  than  the  maximal  overlap  scheme.  The 
observed  difference  nonetheless  exceeds  10  W  m  ‘  over 


Kinorr  t  i  ai,  (,(  m  rcsi’onsi-  lo  c  hanc'.h  in  c  i.oi  d  soi.ar  h)rc  in(, 


a)  Surface 


b)  Atmosphere 


Figure  9.  DitTercnce  (in  watts  per  square  meter)  with  respect  to  the  mean  amount  of  solar  radiation 
absorbed  (a)  at  the  surface  and  (b)  in  the  atmosphere  between  values  obtained  using  the  maximal  overlap 
scheme  and  values  obtained  using  the  grid  cell  averaging  scheme  for  the  NOOAPS  3.3  July  simulations 
I'he  contour  interval  in  Figure  9a  is  20  W  m  ‘  and  that  in  Figure  9b  is  5  W  m  N  The  zero  contour  has  been 
omitted.  Negative  contours  are  indicated  by  dashed  lines.  Shading  is  defined  as  in  Figure  2. 


extensive  regions,  particularly  the  tropical  Pacific  Ocean  and 
the  tropical  and  northern  hemisphere  subtropical  Atlantic 
Ocean. 

5.2.  Response  of  Radiation  Fields 

One  view  of  the  model's  response  to  the  change  in  solar 
radiation  treatments  is  seen  in  Figure  10.  Figure  lOa  shows 
the  difference  between  the  diagnostic  values  from  the 
NOGAPS  3.3  simulations  of  mean  TOA  net  solar  irradiance 
based  on  the  maximal  overlap  treatment  and  the  values 
calculated  in  the  same  simulations  using  the  cloud  grid  cell 
averaging  scheme.  This  plot  represents  the  perturbation  to 
the  total  absorption  (surface  and  atmospheric)  of  solar 
radiation  that  results  from  the  model  change.  In  some 
regions  the  perturbation  is  considerable,  the  peak  difference 
plotted  here  being  104  W  m  '.  For  comparison,  the  differ¬ 
ence  between  the  TO.A  solar  irradiance  field  obtained  with 
NOGAPS  3.3MO  and  the  corresponding  field  obtained  with 
the  cloud  grid  cell  averaging  scheme  in  the  NOGAPS  3.3 
simulations  is  shown  in  Figure  lOb.  The  plots  exhibit  many 
similarities,  both  showing  considerably  greater  net  solar 
irradiance  from  the  maximal  overlap  scheme  for  the  regions 
where  the  simulated  convective  rainfall  rates  are  greatest 


(sec  Figure  3)  There  are  also  some  significant  differences 
betw'een  the  plots,  reflecting  the  effect  of  the  model  response 
to  the  change  in  the  treatment  of  cloud  solar  forcing.  The 
response  over  continental  convective  regions  is  most  nota¬ 
ble.  showing  a  negative  feedback  to  the  perturbation  in  solar 
irradiance  seen  in  Figure  lOa.  This  feedback  results  for  the 
most  part  from  the  increase  in  the  amount  of  convection  over 
land  in  the  NOGAPS  3.3MO  simulations  (see  Figure  3b)  that 
is  generated  by  an  increase  in  the  flux  of  solar  radiation  at  the 
surface.  Such  a  feedback  was  similarly  suggested  in  a  GCM 
experiment  carried  out  by  Rondiill  ct  <il.  [1989],  in  which  the 
effects  of  cumulonimbus  anvils  were  investigated. 

The  model  response  is  more  directly  seen  in  Figures  I  la 
and  lib.  In  Figure  I  la  the  difference  is  plotted  between  the 
values  of  TOA  net  solar  irradiance  predicted  in  the  simula¬ 
tions  with  NOGAPS  3.3MO  and  the  corresponding  results 
diagnosed  from  the  NOGAP.S  3.3  simulations  using  the 
maximal  overlap  scheme.  Figure  I  lb  shows  the  difference  in 
mean  Tf)A  outgoing  longwave  irradianee  between  the  two 
sets  of  simulations.  In  Figure  lla  the  model  response  is 
dominated  by  an  effective  decrease  in  the  amount  of  net 
solar  irradiance  over  continental  regions.  This  decrease 
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a)  NOGAPS  3.3  (maximal  overlap  -  grid-cell-averaging) 


bj  NOGAPS  3.3M0  -  NOGAPS  3,3  (grid-cell-averaging) 


Figure  10.  (a)  Difference  (in  lAuIts  per  square  meter)  between  July  mean  I'O.A  net  solar  irradiance 

diagnosed  from  the  NO(j.-\PS  simulations  using  the  maximal  overlap  solar  radiation  scheme  and  the 
values  obtained  using  the  grid  cell  averaging  scheme  from  the  same  simulations  and  (bl  difference  (in  watts 
per  square  meter)  between  Julv  mean  I'O.A  net  solar  irradiance  from  NOG.APS  3  .^M()  and  the 
corresponding  results  obtained  using  the  grid  cell  tiveraging  treatment  in  the  NOfi.APS  3,3  simulations. 
Shading  and  contours  are  dcHned  as  in  Figure  2. 


reflects  the  negative  feedback  discussed  in  the  preceding 
paragraph,  through  which  the  reduction  in  albedo  that  re- 
ults  from  the  change  in  the  treatment  of  cloud  solar  forcing 
tends  to  generate  an  increased  amount  of  convection  over 
landmasses.  I'he  peak  differences  observed  here,  in  the 
range  of  30-50  W  m  ’  in  some  ureas,  are  fairly  large. 

Whereas  the  solar  irradiance  differences  in  Figure  I  la  are 
more  pronounced  over  continental  regions,  the  correspond¬ 
ing  differences  in  outgoing  longwtive  radiation  seen  in  Figure 
lib  are  more  balanced  in  their  distribution  over  land  ;ind 
water.  There  is  a  moderate  decrease  in  outgoing  longwave 
irradiance  off  the  west  coast  of  .Africa  for  the  simulations 
with  NCXl.APS  3.3MO  and  also  over  a  portion  of  India, 
corresponding  to  similar  solar  irradian -  e  differences  shown 
in  Figure  11a.  These  results  suggest  an  increase  m  the 
amount  of  cloudiness  predicted  for  these  regions  with  the 
revised  model  Increases  in  outgoing  longwave  irradiance 
are  also  observed,  especially  over  the  oceans,  central  .Asia, 
and  the  northwest  coast  of  North  .America,  Whereas  for 
Central  .Americ.i,  much  of  eastern  North  America,  and  parts 
of  the  Pacific  coast  of  .Asia  Figure  11a  shows  some  fairly 
large  differences  in  net  solar  irradiance  associated  with  the 


response  to  the  model  change,  the  differences  in  TO.A 
outgoing  longwave  irradiance  for  these  regions  are  rather 
small. 

The  latitudinal  distribution  of  the  perturbation  in  TO.A 
downwelling  solar  irradiance  (Figure  lOa)  that  rc'iilts  from 
the  change  in  the  treatment  of  cloud  solar  forcing,  and  that  of 
the  model  response  in  terms  of  TO.A  downwelling  sivlar 
(Figure  I  la)  md  longwave  (Figure  lib.  with  a  change  in 
sign)  irradiance.  are  plotted  in  Figure  12.  Results  are  shown 
separately  for  land  and  ice  (Figure  12a)  .md  water  (Figure 
12b).  Although  the  largest  perturbations  occur  over  water, 
the  percentage  response  is  genertilly  greatest  over  land.  In 
the  tropics  over  land,  there  is  a  substantial  negative  response 
in  TO.A  downwelling  solar  irradiance,  tending  to  counter  the 
perturbation,  and  a  positive  response  in  TO.A  downwelling 
longwave  irradiance.  Over  land  in  the  northern  hemisphere 
midlatitudes  the  response  in  terms  of  downwelling  solar 
irnidiance  is  considerably  s.  not  only  in  an  abscilute  sense 
but  also  in  relation  to  the  perturbation.  The  corresponding 
response  in  downwelling  longwave  irradiance  at  midlatitudes 
IS  negative,  in  contrast  to  the  situation  in  the  tropics  The  net 
response  in  TO.A  downwelling  irradiance  in  both  the  iropics 
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a) 


b) 


Hinuri'  1 1.  (ill  DilTcrcncc  (in  wads  per  square  meter)  between  July  mean  I'OA  net  solar  irradiances  from 
the  NOtiAF'S  VJNtO  simulations  and  the  corresponding  results  from  the  NOGAPS  3.3  simulations 
diagnosed  using  the  maxiniid  inerlap  sol;ir  r;idi;ition  scheme  and  (b)  difference  /in  walls  per  square  rneler) 
heiueen  .lulv  mean  fOA  outgoing  long\c;ive  irradiances  obtained  in  the  simulations  with  NOGAPS  3.3MO 
and  NOGAPS  3  3  The  cimtour  intersiil  in  both  figures  is  10  W  m  ■,  and  the  zero  contour  has  been 
omitied.  Negative  contours  ;irc  indicated  by  diished  lines.  Shading  o  defined  as  in  Figure  2. 


<md  midlatitudes  is  negative  Despite  the  significantly  greater 
response  in  fO  A  dounwellmg  solar  irradiance  over  land  in 
the  tropic  the  net  response  in  fO.A  downwelling  irradiance 
as  a  percentage  of  the  mean  perturbation  tends  to  he  greater 
in  magnitude  at  midlatitiides.  as  seen  in  fable  I. 

6.  Response  of  Model  Cloud  Fields 

fhe  radiation  differences  seen  in  Figure  1 1  result  largely 
from  changes  in  the  predicted  cloud  fields.  I'hese  changes 
m.is  include  changes  both  in  cloud  fraction  and  in  cloud 
optical  thickness  at  different  model  grid  levels.  The  change 
in  total  cloiiiliness  that  results  from  the  change  in  the 
treatment  of  cloud  solar  forcing  is  shown  in  Figure  13.  The 
increases  in  cloudiness  off  the  west  coast  of  Africa  and  over 
a  portion  of  India  that  are  suggested  by  the  change  in 
outgoing  longwave  radiation  in  Figure  lib  are  both  evident 
in  this  plot.  In  general,  the  changes  in  total  cloudiness  tend 
to  be  better  correhited  with  changes  in  longwave  rather  than 
solar  irradiance.  .Although  the  changes  in  cloudiness  ob¬ 
served  over  land  are  generally  associated  with  changes  in 
convective  activity,  the  changes  over  the  oceans  tend  to 


occur  away  from  convective  regions  (see  Figure  3)  and 
generally  reflect  a  decrease  in  moisture  at  mid  to  uppe'' 
levels.  1  he  significant  decreases  in  cloudiness  at  midlati¬ 
tudes  are  responsible  for  the  negative  response  of  the  model 
there  in  terms  of  TOA  longwave  irradiance  noted  in  section 
5.2. 

•Although  cloud  optical  thickness  data  were  not  saved  from 
the  model  runs,  estimates  cif  cloud  optical  thickness  were 
made  on  the  basis  of  model  July  mean  temperature  fields  for 
cloud  layers  .^(M)  m  thick.  The  fractional  changes  in  optical 
thickness  obtained  for  the  two  model  versions  at  8.^0  mbar. 
.‘'00  mbar,  and  2(K)  mbar  are  plotted  in  Figure  14.  We  find  that 
cloud  optical  thicknesses  are  generally  greater  for  NGGAPS 
3,3MO.  reflecting  a  warming  of  the  troposphere  that  w  ill  be 
discussed  in  section  7.3.  Note  that  only  results  obtained  for 
stable  clouds  are  shown  here.  Results  obtained  for  convec¬ 
tive  clouds  are  generally  similar,  except  w  here  temperatures 
are  substantially  below  freezing.  Wheie  clouds  composed  of 
both  liquid  water  and  ice  exist,  the  fractional  changes 
obtained  for  convective  clouds  are  generally  considerably 
smaller  than  those  for  stable  clouds  because  of  the  greater 
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Figure  12.  Zonal  perturbation  in  TOA  downwelling  solar 
irradiance  represented  b\  the  change  in  the  treatment  of 
cloud  sohir  forcing  (plus  signs l  and  changes  in  TOA  down- 
welling  solar  (circles)  and  longwave  (asterisks)  irradiance 
that  result  from  the  model  response  to  that  perturbation. 
Results  are  shown  R'r  (a)  land  and  ice  ;ind  (b)  water. 


ice  contents  assumed  for  convective  clouds.  At  lower  levels 
the  increases  in  optical  thickness  occur  largely  over  conti¬ 
nental  regions,  though  generally  not  over  the  most  convec- 
tively  active  areas,  as  seen  by  comparison  with  Figure  .^.  The 
increases  in  cloud  optical  thickness  tend  to  counter  the 
increased  surface  absorption  of  solar  radiation  that  results 
from  the  change  in  the  treatment  of  cloud  solar  forcing. 


Table  I.  C  hanges  in  fOA  Downwelling  Irradiance  Over 
Northern  Hemisphere  l.and 


\  .atitiKlmal 
/(inc 

Pertiirhi^lion. 

W  ni  * 

Rc 

Solar, 

W  m 

bptinsc 

longwave. 
W  m  ■ 

I'otal. 

'  } 

0'  -9  N 

4*;  s 

10  1 

(  )  .40.0 

9-IS  N 

tx  -! 

14  9 

X  0 

(  )  18,0 

IX  -27  N 

10. X 

X.2 

1  )  ll.s 

27  -76  N 

IX  t 

X..S 

17 

(  I  .47. .S 

t6  -IS  N 

21)  (1 

^  7 

VI 

(  )  44, X 

4S  -S4  \ 

2X  t 

s  X 

7  4 

(  )  46.7 

S4 -6t  N 

X  4 

4  0 

(  )  .48.7 

6t  -72  N 

19  ^ 

2.4 

(  )  29  1 

72  -XI  N 

|S  t 

It) 

1  s 

1  )  16.5 

XI  -WN 

9  6 

17 

4  0 

24.7 

though  their  actual  effect  on  the  solar  radiation  field  will 
depend  on  the  vertical  profile  of  cloudiness  in  a  given  region. 

The  effects  of  changes  in  cloudiness  on  the  radiation  fields 
are  most  clearly  seen  in  the  corresponding  changes  in  cloud 
radiative  forcing.  Plots  showing  the  changes  in  surface  cloud 
solar  forcing  and  atmospheric  cloud  longwave  forcing  that 
result  from  incorporating  the  maximal  overlap  scheme  into 
the  model  are  shown  in  Figures  l.'Sa  and  l.‘'b,  respectively. 
The  cloud  solar  forcing  changes  shown  in  Figure  l.'ia  were 
calculated  using  the  maximal  overlap  treatment  for  both  sets 
of  simulations.  Negative  values  in  Figure  15a  reflect  in¬ 
creased  cloud  solar  forcing,  and  positive  values  in  Figure  15b 
reflect  increased  cloud  longw'ave  forcing.  We  do  not  show 
the  changes  in  cloud  solar  forcing  of  the  atmosphere  or  cloud 
longwave  forcing  of  the  surface,  because  these  changes  are 
considerably  smaller  than  those  plotted  here,  as  one  may 
expect  on  the  basis  of  the  modeling  study  of  Slinf;o  and 
Slir.fio  (1988).  To  a  considerable  degree,  changes  in  surface 
cloud  solar  forcing  are  correlated  with  perturbations  in 
surface  solar  absorption  shown  in  Figure  9a.  Although  the 
largest  such  perturbations  occur  in  the  tropics,  they  occur 
most  extensively  there  over  the  oceans  where  there  is  no 
surface  response  in  the  present  model;  hence  the  observed 
cloud  feedback  over  land  is  more  evenly  distributed  latitu- 
dinally  than  might  otherwise  be  expected.  The  relative 
response  of  the  model  in  terms  of  the  ratio  of  the  change  in 
cloud  surface  forcing  to  the  perturbation  shown  in  Figure  9a 
is  in  many  regions  very  large,  showing  the  importance  of 
improvements  in  the  treatment  of  cloud  solar  forcing  for 
cloud  field  prediction.  Over  the  oceans  the  effect  of  cloud 
feedback  on  absorption  of  longwave  radiation  by  the  atmo¬ 
sphere  tends  to  reinforce  the  direct  cooling  effect  of  the 
change  in  the  treatment  of  cloud  solar  forcing  seen  in  Figure 
9b.  These  effects  are  not  strongly  correlated  spatially,  how¬ 
ever,  there  being  a  greater  effect  of  cloud  feedback  at  higher 
latitudes  than  might  be  expected  from  the  direct  perturbation 
in  atmospheric  absorption.  Over  land,  because  of  the  in¬ 
creased  convective  activity  that  results  from  the  surface 
response,  changes  in  atmospheric  absorption  of  longwave 
radiation  tend  to  counter  the  direct  cooling  effect  seen  in 
Figure  9b.  The  latitudinal  dependencies  of  the  changes  in 
cloud  solar  forcing  of  the  surface  and  longwave  forcing  of  the 
atmosphere  are  shown  in  Figure  16  for  land  and  ice  (Figure 
16a)  and  water  (Figure  16b).  A  comparison  with  Figure  12 
shows  that  most  of  the  change  in  TO.A  downwelling  solar 
irradiance  shown  there  is  accounted  for  by  the  changes  in 
surface  absorption  shown  in  Figure  16.  To  some  extent, 
changes  in  TOA  downwelling  longwave  irradiance  and 
changes  in  atmospheric  absorption  are  similarly  related,  hut 
in  this  case,  changes  in  temperature  and  moisture  also  play  a 
significant  role. 

7.  Response  of  Model  Precipitation 
and  Temperature  Fields 

7.1.  Precipitation  Response 

The  change  in  convective  precipitation  rate  was  shown  in 
Figure  .56.  Although  most  of  the  precipitation  increases  over 
land  arc  from  convective  precipitation,  there  is  a  large 
increase  (as  much  as  8  mm  d  ')  in  large-scale  precipitation 
over  land  near  the  Bay  of  Bengal.  The  change  in  the 
difference  between  precipitation  and  evaporation  is  plotted 
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Figure  13.  Change  in  July  mean  total  cloudiness  (assuming  maximally  overlapping  clouds).  The  contour 
interval  is  0.05.  Shading  is  defined  as  in  Figure  2. 


in  Figure  17.  showing  that  there  is  a  tendency  toward  an 
increase  in  this  quantity  over  several  continental  regions, 
most  notably  near  the  Bay  of  Bengal  and  central  Africa, 
reflecting  an  increased  flux  of  moisture  from  sea  to  land.  The 
relationship  between  convective  rainfall  rate  and  surface 
cloud  solar  forcing  is  also  affected  by  the  model  change.  This 
fact  is  illustrated  in  Figure  18.  which  shows  for  the  two 
model  versions  the  surface  cloud  solar  forcing  plotted 
against  convective  precipitation  rate  for  9°  x  9°  continental 
regions  extending  from  0'  to  9‘’N.  We  find,  as  one  expects, 
that  convective  rainfall  tends  to  increase  with  increased 
surface  cloud  solar  forcing  and  that  for  a  given  convective 
rainfall  rate,  the  surface  cloud  solar  forcing  predicted  by 
NOGAPS  generally  exceeds  that  predicted  by  NOGAPS 
.V.fMO  by  .70-50  W  m  v  Somewhat  greater  changes  in  this 
relationship  are  observed  over  water.  These  changes  are 
reflected  over  land  in  changes  in  the  corresponding  relation¬ 
ships  for  other  surface  fluxes. 

7.2.  Surface  Response 

As  shown  in  Figure  19.  there  are  some  significant  ground 
temperature  changes  associated  with  the  change  in  the 
treatment  of  cloud  solar  forcing.  Ground  temperatures  gen¬ 
erally  increase  over  continental  regions,  consistent  with 
what  one  may  expect  on  the  basis  of  the  increase  in  surface 
solar  irradiance  shown  in  Figure  9a.  The  relationship  be¬ 
tween  the  temperature  increase  and  the  perturbation  in  the 
absorption  of  solar  radiation  at  the  surface  shown  in  Figure 
9a,  however,  is  found  to  vary  considerably.  Note,  in  partic¬ 
ular,  that  the  large  increase  in  ground  temperature  over 
central  Asia  occurs  where  the  perturbation  in  surface  solar 
irradiance  is  less  than  40  W  m  ’.  whereas  larger  perturba¬ 
tions  over  central  Africa,  India,  and  North  and  South 
America  are  accompanied  by  significantly  smaller  tempera¬ 
ture  increases.  The  observed  response  of  the  hind  surface 
must  therefore  reflect  variations  in  the  mode!  response  in 
terms  of  cloud  feedback  and  latent  and  sensible  heating.  One 
finds,  for  example,  by  comparing  these  observations  with 
Figure  15a,  where  the  effect  of  cloud  feedback  on  surface 
solar  irradiance  is  shown,  that  the  region  over  central  Asia 
where  there  is  a  substantial  increase  in  ground  temperature 
is  distinguished  by  the  absence  of  the  negative  feedback 
observed  over  other  regions.  There  appears  to  be  some 


tendency,  in  fact,  for  the  model  to  produce  a  positive 
feedback  in  this  region.  Referring  to  Figure  3b.  one  sees  that 
the  general  increase  in  convective  precipitation  over  conti¬ 
nental  regions  is  not  observed  over  much  of  central  .Asia, 
possibly  the  result  of  a  change  in  the  large-scale  circulation. 

Although  the  lack  of  a  negative  cloud  feedback  on  surface 
solar  irradiance  over  central  Asia  helps  to  explain  the 
variability  in  the  simulated  changes  in  ground  temperature  to 
some  extent,  this  facet  of  the  simulations  does  not  fully 
explain  the  simulated  temperature  change  variations.  This 
fact  can  be  seen  by  comparing  the  ground  temperature 
differences  in  Figure  19  with  Figime  20  which  shows  the  net 
change  in  surface  solar  irradiance.  including  the  effect  of  the 
model  response.  The  net  change  in  surface  solar  irradiance  is 
to  a  considerable  degree  a  worse  indicator  of  ground  tem¬ 
perature  changes  than  the  perturbation  in  surface  solar 
irradiance  plotted  in  Figure  9a.  For  example,  although  the 
net  increase  in  surface  solar  irradiance  over  portions  of  India 
exceeds  the  corresponding  increase  over  central  .Asia,  the 
change  in  ground  temperature  over  India  is  very  small.  On 
the  other  hand,  the  observed  increase  in  ground  temperature 
over  central  Africa  is  accompanied  by  only  scattered  modest 
increases  in  surface  solar  irradiance. 

The  lack  of  any  significant  increase  in  ground  temperature 
over  India  is  explained  by  Figure  21.  which  shows  the 
change  in  surface  latent  heating  (positive  values  indicating 
increased  evaporation  and  surface  cooling).  There  is  a  region 
of  significantly  enhanced  evaporation  over  India  which  over¬ 
laps,  but  is  not  coincident  with,  a  region  of  enhanced 
convective  precipitation  seen  in  Figure  3b.  This  increase  in 
surface  evaporation  largely  offsets  the  increase  in  surface 
solar  irradiance  over  India  seen  in  Figure  20.  The  change  in 
surface  evaporation  is  almost  everywhere  considerably 
greater  than  either  the  change  in  surface  sensible  heating  or 
the  change  in  the  net  surface  emission  of  longwave  radiation. 
F,xceptions  occur,  however,  in  arid  regions  such  as  central 
Asia,  where  increases  in  evaporation  are  necessarily  limited. 
The  increased  evaporation  appears  to  result,  to  a  consider¬ 
able  extent,  from  increased  surface  wind  velocities,  though 
in  regions  of  enhanced  precipitation,  increased  low-level 
drying  caused  by  the  model's  Arakawa-.Schubert  cumulus 
parameterization  and  increased  availability  of  surface  mois- 
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a) 


b) 


c) 


Figure  14.  Fractional  change  in  optical  thickness  of  stable  clouds  calculated  using  July  mean  tempera¬ 
ture  fields.  Results  are  shown  for  (a)  850  mbar.  tb)  500  mbar,  and  (c)  200  mbar.  The  contour  interval  is  0. 1 . 
Shading  is  defined  as  in  Figure  2. 


ture  for  evaporation  are  also  likely  factors.  Over  the  oceans, 
there  appears  to  be  a  tendency  for  increased  evaporation 
where  there  has  been  a  reduction  in  low-level  moisture. 
Comparing  Figure  21  with  Figure  22,  which  shows  the 
change  in  1000-mbar  specific  humidity,  we  see  such  a 
correlation,  for  example,  to  the  east  of  Madagascar  and  off 
the  west  coast  of  South  America, 


7.3.  Atmospheric  Stability 

The  observed  increases  in  convective  activity  reflect  de- 
stabilization  of  the  atmosphere  near  the  surface  brought 
about  by  the  increased  surface  heating  and  the  response  of 
the  model's  revised  Arakawa-Schubert  convective  parame- 
terizdtion.  Changes  in  the  temperature  lapse  rate  of  the 
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a) 


b) 


Figure  15.  Change  in  July  mean  surface  cloud  (a)  solar  forcing  and  atmospheric  cloud  tb)  longwave 
forcing.  The  results  plotted  in  Figure  15a  wjere  obtained  using  the  maximal  overlap  treatment  for  both  sets 
of  simulations.  The  contour  interval  in  both  plots  is  5  W  m”’.  Shading  is  defined  as  in  Figure  2. 


atmosphere  over  convective  regions  show  to  some  extent 
the  degree  to  which  the  convective  parameterization 
counters  the  tendency  toward  decreased  atmospheric  stabil¬ 
ity.  though  changes  in  atmospheric  moisture  can  signifi¬ 
cantly  affect  atmospheric  stability  as  well.  Concerning  an 
experiment  w  ith  the  UCL.A/GLA  GCM.  Randall  '•t  al.  f  1989) 
hypothesized  that  the  cause  of  a  substantial  increase  in 
temperature  in  the  northern  hemisphere  midlatitudes  that 
tended  to  stabilize  the  atmosphere  was  the  response  of  the 
Arakawa-Schubert  sch'^me  in  the  model  to  destabilization 
brought  about  by  increased  atmospheric  moisture.  In  the 
Randall  et  al.  experiment  the  GCM  was  run  for  July  condi¬ 
tions,  and  the  effect  of  convective  anvils  was  tested  by 
removing  them.  The  model  change  studied  in  the  present 
work  has  a  similar,  though  not  as  dramatic,  effect  on  solar 
radiation  fields  as  removing  convective  anvils.  There  are 
some  significant  differences  between  these  studies  as  well, 
however,  including  the  fact  that  in  the  present  work  there  is 
no  direct  effect  on  longwave  radiation,  and  cloud  feedback  is 
allowed.  In  order  to  investigate  any  changes  in  atmospheric 
conditional  instability  associated  with  changes  in  tempera¬ 
ture  in  the  present  simulations,  we  show  in  Figure  23  the 
change  in  temperature  lapse  rate  betweeii  the  850-mbar  and 
the  500-mbar  levels.  We  see  some  evidence  of  stabilization 
(decreased  lapse  rate)  at  mid  to  high  latitudes  in  the  northern 
hemisphere,  though  there  are  regions  of  destabilization  in 
these  latitudinal  zones  as  well,  particularly  over  Asia  and 
central  North  America.  Comparing  this  plot  with  the  plot  of 


changes  in  the  convective  precipitation  rate  shown  in  Figure 
3b.  we  note  that  the  increase  in  convective  precipitation  over 
Colombia  seen  in  Figure  3b  corresponds  to  a  small  increase 
in  temperature  lapse  rate,  whereas  increased  convective 
precipitation  over  western  Aifrica  occurs  where  there  is  a 
small  decrease  in  temperature  lapse  rate.  Differences  in 
atmospheric  moisture  may  explain  this  difference,  but  this 
point  is  not  clear.  A  look  at  Figure  24.  which  shows  the 
change  in  precipitable  water  amount,  reveals  a  significant 
increase  over  western  Africa,  consistent  with  the  theory  that 
destabilization  by  increased  low-level  moisture  in  convec¬ 
tive  regions  results  in  a  decrease  in  temperature  lapse  rate. 
Over  Colombia  there  is  a  large  region  with  no  significant 
change  in  precipitable  water  (suggesting  that  the  increases  in 
temperature  lapse  rate  in  that  region  is  not  inconsistent  with 
this  theory),  but  a  regional  increase  in  precipitable  water  is 
also  observed  there.  It  should  be  noted  that  some  of  the 
largest  changes  in  lapse  rate  in  Figure  23  actually  occur  aw  ay 
from  convectively  active  regions,  supporting  the  argument 
against  the  model's  convective  scheme's  being  the  primary 
cause  of  such  changes. 

The  extent  of  atmospheric  temperature  changes  is  more 
clearly  seen  in  Figure  25.  which  shows  the  temperature 
changes  at  850  mbar  (Figure  25a).  500  mbar  (Figure  25b).  and 
200  mbar  (Figure  25c).  The  large  increase  in  temperature  at 
midlatitudes  in  the  northern  hemisphere  bears  some  resem¬ 
blance  to  the  zonal  temperature  increase  noted  by  Randall  cl 
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Figure  16.  Zonal  changes  in  surface  absorption  of  solar 
radiation  (circles)  and  atmospheric  absorption  of  longwave 
radiation  (asterisks)  that  result  from  the  model  response  to 
the  change  in  the  treatment  of  cloud  solar  forcing.  Results 
are  shown  for  (a)  land  and  ice  and  (b)  water. 


III.  [1989],  The  temperature  increase  at  850  mbar  observed 
here  is  greatest  over  arid  central  Asia,  where  the  large 
ground  temperature  increase  noted  in  the  preceding  section 
is  observed  (Figure  19).  There  is  also  a  large  region  of 
increased  temperature  over  northern  North  America.  At  500 
mbar  the  region  of  increased  temperature  is  considerably 
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Figure  18.  July  mean  surface  cloud  solar  forcing  as  a 
function  of  convective  precipitation  rate  for  9°  x  9°  conti¬ 
nental  regions  extending  from  0°  to  9°N.  Results  are  shown 
for  NOGAPS  3.3  (plus  signs)  and  NOGAPS  3.3MO  (circles). 
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more  extensive  but  largely  confined  to  mid  and  high  lati¬ 
tudes.  At  200  mbar  the  greatest  temperature  increases  are 
found  over  the  Atlantic  and  northern  Pacific  Oceans.  Noting 
the  strong  correlation  between  increases  in  850-mbar  tem¬ 
perature  and  increases  in  ground  temperature,  it  appears  that 
the  low-level  temperature  increases  in  the  present  simula¬ 
tions  result  from  increases  in  surface  sensible  heating  and 
emission  of  longwave  radiation  and  from  increases  in  verti¬ 
cal  eddy  or  shallow  convective  transport.  The  temperature 
increases  at  mid  to  upper  levels  may,  to  some  extent,  reflect 
the  vertical  eddy  transport  and  advection  of  air  heated  at 
lower  levels.  Changes  in  radiative  heating  may  possibly  be  a 
factor  in  the  temperature  changes  as  well.  Although  the 
direct  effect  of  the  model  change  on  atmospheric  heating, 
shown  in  Figure  9b,  is  to  cool  the  atmosphere.  Figure  15b 
shows  that  the  cloud  response  to  the  model  change  tends  to 
warm  the  atmosphere  over  certain  regions,  including  south¬ 
ern  Asia  and  western  Africa. 

8.  Discussion  and  Conclusions 

The  NOGAPS  forecast  model  results  and  ERB  satellite 
data  presented  in  section  4  show  a  considerable  improve¬ 
ment  in  solar  radiation  fields  with  the  assumption  of  maxi- 


Figure  17.  Change  in  July  mean  difference  between  the  rate  of  precipitation  and  the  rate  of  evaporation. 
The  contour  interval  is  2  mm  d  ' .  Shading  is  defined  as  in  Figure  2. 
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Figure  19.  Change  in  July  mean  surface  temperature  (in  degrees  Kelvin).  The  contour  interval  is  I  K. 
Shading  is  defined  as  in  Figure  2. 


mally  overlapping  clouds  in  the  treatment  of  cloud  solar 
forcing.  The  improved  correlation  between  the  TOA  solar 
and  the  longwave  irradiances  demonstrated  in  section  4.3 
supports  the  conclusion  that  this  improvement  is  not  merely 
the  result  of  a  compensation  for  errors  due  to  inadequacies  in 
the  model's  cloud  parameterization.  The  results  obtained 


should  be  expected,  because  the  change  in  the  solar  radia¬ 
tion  scheme  addressed  here  brings  consistency  between  the 
treatment  of  fractional  cloudiness  in  the  solar  and  longwave 
radiation  schemes  of  the  model.  Because  of  the  relatively 
low  optical  thickness  of  most  nonconvective  clouds  param¬ 
eterized  in  the  GCM  simulations  studied  here,  the  most 
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Figure  20.  Change  in  July  mean  surface  solar  irradiance  (in  watts  per  square  meter),  reflecting  both  the 
change  in  cloud  solar  forcing  treatment  and  the  model  response.  The  contour  interval  is  20  W  m  -. 
Shading  is  defined  as  in  Figure  2. 
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Figure  21.  Change  in  July  mean  surface  latent  heating  (in  watts  per  square  meter).  The  contour  interval 
is  10  W  m  Shading  is  defined  as  in  Figure  2. 
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Figure  22.  Change  in  July  mean  lOOO-mbar  specific  humidity  (in  grams  per  kilogram).  The  contour 
interval  is  0.2  g  kg~'.  Shading  is  defined  as  in  Figure  2. 

significant  improvement  in  model  solar  radiation  fields  was  over  continental  regions  that  results  from  the  increased 

shown  in  section  4  to  be  in  regions  of  deep  convection.  We  amount  of  surface  solar  radiation.  Along  with  the  increase  in 

found,  for  example,  that  improvements  in  simulated  surface  continental  convection,  there  is  a  tendency  for  an  increased 

solar  irradiances  exceed  100  V.'  m"’  in  portions  of  the  flux  of  moisture  from  sea  to  land.  Substantial  fractional  in- 

tropical  Pacific  Ocean.  creases  in  cloud  water  content  are  found  to  occur,  for  the  most 

The  dominant  model  response  in  the  present  study  is  the  part  associated  with  extensive  regions  of  increased  tropo- 
observed  increase  in  convective  precipitation  and  cloudiness  spheric  temperature  in  the  northern  hemisphere  midlatitudes. 
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Figure  23.  Change  in  July  mean  temperature  lapse  rate  between  850-mbar  and  500-mbar  levels.  The 
contour  interval  is  0.1  K  km  '.  Shading  is  defined  as  in  Figure  2. 


Figure  24.  Change  in  July  mean  precipitable  water  amount  (in  millimeters).  The  contour  interval  is  2 
mm.  Shading  is  defined  as  in  Figure  2. 
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a) 


Figure  25.  Change  in  July  mean  temperature  (in  degrees  Kelvin)  at  (a)  850  mbar,  (b)  500  mbar.  and  (c) 
200  mbar.  The  contour  interval  is  I  K.  Shading  is  defined  as  in  Figure  2. 


Clouds  generally  act  as  a  negative  feedback  to  the  pertur¬ 
bation  in  solar  irradiance  represented  by  the  model  change. 
Estimates  of  this  perturbation  were  obtained  as  the  differ¬ 
ence  between  solar  irradiances  diagnosed  (but  not  used  in 
the  model)  assuming  maximally  overlapping  clouds  and  solar 
irradiances  produced  (and  used  in  the  model)  in  the  same 
simulations  by  the  unrevised  NOGAPS  cloud  solar  forcing 
scheme  using  grid  cell  averaged  clouds.  With  respect  to  the 
perturbation  defined  in  this  way,  the  model  response  in 


terms  of  the  change  in  surface  cloud  solar  forcing  computed 
for  maximally  overlapping  clouds  is  in  some  regions  quite 
large.  This  large  response  is  responsible  for  the  relatively 
small  reduction  observed  in  the  error  in  TOA  downwelling 
solar  irradiance  over  land  in  certain  regions,  notably  India 
and  portions  of  Africa  and  eastern  North  America,  and 
shows  the  importance  of  improvements  in  the  treatment  of 
cloud  solar  forcing  for  cloud  field  prediction. 

Over  land  in  the  northern  hemisphere  the  model  response 
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in  terms  of  the  change  in  TOA  downwelling  solar  irradiance 
evaluated  using  the  maximal  overlap  treatment  is  apprecia¬ 
bly  greater  in  the  tropics  than  at  midlatitudes.  For  example, 
between  the  equator  and  9'’N  the  mean  response  is  -23.7  W 
m  or  52‘'r  of  the  45.4  W  m  “  perturbation  in  TOA 
downwelling  solar  irradiance  represented  by  the  model 
change.  The  corresponding  response  between  45°N  and  54°N 
is  -5.8  W  m  -.or  215?  of  the  28,3  W  m  "  perturbation  at 
that  latitudinal  zone.  We  find,  however,  that  because  of 
latitudinal  difl'erences  in  the  response  of  the  model's  long¬ 
wave  radiation  field  the  percentage  net  response  in  TOA 
downwelling  irradiance  at  midlatitudes  tends  to  exceed  that 
in  the  tropics. 

Changes  in  ground  temperature  observed  here  reflect  a 
number  of  factors,  including  cloud  cover,  surface  moisture, 
and  large-scale  dynamics.  The  tendency  for  ground  temper¬ 
ature  to  increase  in  convective  regions  because  of  increased 
surface  solar  irradiance  is  moderated  to  a  considerable 
degree  by  cloud  feedback  and  increased  surface  evapora¬ 
tion.  Increases  in  surface  sensible  heating  and  net  emission 
of  longwave  radiation  are  limited  largely  to  regions  where 
these  other  factors  do  not  act  to  limit  temperature  increases. 
Central  Asia  is  the  most  notable  example  of  such  regions 
observed  in  the  present  study.  It  experiences  a  decrease  in 
precipitation  and  cloud  cover.  Evaporation  does  not  in¬ 
crease  enough  to  be  a  significant  factor  in  limiting  the 
resulting  ground  temperature  increase  because  of  the  small 
amount  of  ground  moisture.  The  extensive  atmospheric 
temperature  increases  in  the  northern  hemisphere  midlati¬ 
tudes  appear  to  be  associated,  at  least  at  low  levels,  with 
such  ground  temperature  increases. 

.Although  imp'sived  agreement  of  NOGAPS  solar  radiation 
fields  with  ERB  data  has  been  demonstrated  with  the  max- 
imtil  overtiip  treatment  of  cloud  solar  forcing,  considerable 
errors  remain.  As  indicated  in  section  1.  these  errors  may 
derive  from  a  number  of  sources.  We  have  showm  in  the 
present  study  that  the  model's  solar  radiation  field  is  very 
sensitive  to  the  tretitment  of  fractional  cloudiness,  and 
further  work  is  needed  to  determine  the  extent  of  the 
remaining  errors  associated  with  this  treatment.  On  the  basis 
of  their  study  of  cloud  overlap  statistics  over  the  North 
Atlantic  during  January  1979.  Tian  and  Cnrn-  (1989)  noted 
that  although  the  maximal  overlap  assumption  appears  to  be 
a  good  iipproximation  for  adjacent  cloudy  layers,  the  random 
overliip  assumption  is  a  better  approximation  for  cloudy 
layers  separated  by  clear  interstices.  The  solar  radiation 
scheme  of  (Icicvn  and  Hallini’sworth  (1979)  is  somewhat 
consistent  with  these  observations  (cf.  Ritter  and  Cieleyn. 
I992|.  though  its  tipplication  of  the  maximal  overlap  assump¬ 
tion  for  adjacent  cloudy  layers  assumes  no  horizontal  vari¬ 
ation  of  solar  irraditince  within  the  cloudy  or  clear  portions 
of  ;i  given  model  grid  cell.  Composite  schemes  of  this  type 
ma>  prove  to  be  acceptable  solutions  for  the  radiation 
treatment  of  fractional  cloudiness  in  global  models,  but  a 
compelling  empirical  basis  for  ;my  such  solution  must  await 
further  studies  of  the  regional  and  seasonal  dependencies  of 
the  structure  of  cloud  fields. 

Appendix:  Parameterization  of  Cloud 
Optical  Thickness  in  NOGAPS 

Cloud  optical  thicknesses  produced  by  NOG.APS  3.3  are 
based  on  separate  calculations  for  stable  and  convective 


clouds.  In  general,  the  optical  thickness  is  taken  to  be  a 
weighted  sum  of  values  calculated  for  water  and  ice  clouds. 
The  weighting  varies  linearly  with  temperature  within  the 
temperature  range  233.16-273.16  K  and  produces  optical 
thicknesses  consistent  with  all  ice  clouds  at  temperatures  of 

233.16  K  and  less  and  all  water  clouds  at  temperatures  of 

273.16  K  and  greater. 

For  water  clouds  the  optical  thickness  t/  of  a  cloud  layer 
corresponding  to  a  model  grid  level  of  thickness  Ac  is  taken 
to  be  proportional  to  the  liquid  water  path  following 
Stephens  [1978al.  with  liquid  water  content  vv,  assumed  to 
be  proportional  to  the  saturation  specific  humidity  t/,.  Thus 

T,=  3LWT/2r,,.  (1) 

where 

LWP=wi\z  (2) 

and 

If;  =  ft/,.  (3) 

Appropriate  values  for  the  equivalent  radius  r^  of  the  drop 
size  distribution  and  for  the  parameter  c  vary  widely, 
depending  on  cloud  type.  In  the  current  NOGAPS  parame¬ 
terization  a  value  of  0.0079  is  used  for  c  for  stable  clouds, 
which  yields  a  liquid  water  content  of  0.14  g  m  at  20x'. 
This  liquid  water  content  corresnond.s  to  one  of  the  stratoc- 
umulus  drop  size  distnbutions  in  Table  4  of  Stingo  and 
Sch.eeker  (1982).  The  equivalent  radius  given  for  that  drop 
size  distribution  is  5.4  /xm.  which  is  the  value  used  in 
NOGAPS  forr,,  in  the  case  of  stable  clouds.  For  convective 
water  clouds  the  value  used  for  e  is  0.013.  which  gives  a 
liquid  water  content  of  0.23  g  m  at  20°C.  The  equivalent 
radius  for  this  case  was  taken  to  be  5.5  ^m,  the  same  as  that 
used  by  Stephens  [1978b]  for  fair  weather  cumulus.  The 
liquid  water  content  cited  here  is  less  than  the  value  used  by 
Stephens  (I.O  g  m"’)  and  appears  to  be  too  low,  except  for 
very  small  clouds.  Warner  (1955.  I969J  frequently  observed 
peak  water  contents  in  cumulus  clouds  in  the  range  of 
0. 5-1.5  g  m  '.  though  values  on  the  order  of  259F  less  than 
the  peak  values  may  be  justified  for  use  in  a  cloud  parame¬ 
terization  on  the  basis  of  the  substantial  variability  he 
observed  within  clouds.  Warner  (1955)  generally  found  the 
water  content  to  peak  near  cloud  top,  a  feature  which  is  not 
reproduced  in  the  present  formulation  of  the  scheme  because 
of  the  unrealistic  decline  in  liquid  water  content  with  height 
given  by  (3). 

For  ice  clouds  the  optical  thickness  r,  of  a  given  cloud 
layer  corresponding  to  a  model  grid  level  of  thickness  Ac  is 
obtained  from  its  definition, 

T;  =  Wjkjl:..  (41 

where  >< ,  is  the  ice  water  content  and  A,  is  the  extinction 
coefficient,  fhe  extinction  coefficient  is  given  by 

k,  =  tr,JViPi.  (5) 

where  o-,.  is  the  extinction  cross  section.  V,  is  the  volume  of 
a  typical  ice  crystal,  and  p,  is  the  density  of  ice.  The  value 
used  for  V,  corresponds  to  columnar  ice  crystals  of  radius  60 
Xi,m  and  length  300  /Ltm.  This  crystal  size  falls  within  the 
range  observed  in  several  observational  studies  [Lion.  1986), 
The  corresponding  extinction  cross  section,  5,653  x  10  * 
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m",  given  by  Lion  [1986]  and  based  on  a  numerical  calcula¬ 
tion,  is  used  for  a,,.  Values  used  for  the  ice  water  conte  it 
differ  for  stable  clouds  and  convective  clouds.  For  stable 
clouds,  w ,  is  given  in  units  of  grams  per  cubic  meter  by 

w,  =  0.028  :r|s20°C 

In  (if,)  =  -7.6  +  4  exp  [-0.2443  x  10  ’(jFj  -  20.)-  ']  (6) 

iri  >  20“c 

where  T  is  the  temperature  in  degrees  Celsius.  The  expres¬ 
sion  used  for  |7  greater  than  20°C  approximates  the  expres¬ 
sion  given  by  Liou  [I986|.  which  is  based  on  data  presented 
by  Heymsfield  and  Plait  [19841.  For  |7'1  less  than  20°C  the 
constant  value  assumed  for  w,  gives  an  approximately 
continuous  temperature  dependence  of  w;.  For  convective 
clouds.  If,  is  assigned  a  constant  value  of  0.1  g  m“^,  which 
may  be  somewhat  low.  on  the  basis  of  the  range  of  from  0.15 
to  0.25  g  m  ’  observed  by  Heymsfield  and  Knollenberg 
(1972)  in  cirrus-generating  cells. 
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